Dose and DNA damage modelling of diffusing alpha-emitters radiation therapy using Geant4.

PURPOSE: Diffusing alpha-emitters radiation therapy (DaRT) is a brachytherapy technique using α-particles to treat solid tumours. The high linear energy transfer (LET) and short range of α-particles make them good candidates for the targeted treatment of cancer. Treatment planning of DaRT requires a good understanding of the dose from α-particles and the other particles released in the 224Ra decay chain. METHODS: The Geant4 Monte Carlo toolkit has been used to simulate a DaRT seed to better understand the dose contribution from all particles and simulate the DNA damage due to this treatment. RESULTS: Close to the seed α-particles deliver the majority of dose, however at radial distances greater than 4 mm, the contribution of ß-particles is greater. The RBE has been estimated as a function of number of double strand breaks (DSBs) and complex DSBs. A maximum seed spacing of 5.5 mm and 6.5 mm was found to deliver at least 20 Gy RBE weighted dose between the seeds for RBEDSB and RBEcDSB respectively. CONCLUSIONS: The DNA damage changes with radial distance from the seed and has been found to become less complex with distance, which is potentially easier for the cell to repair. Close to the seed α-particles contribute the majority of dose, however the contribution from other particles cannot be neglected and may influence the choice of seed spacing.

Asymptotic analysis of particle cluster formation in the presence of anchoring sites.

The aggregation or clustering of proteins and other macromolecules plays an important role in the formation of large-scale molecular assemblies within cell membranes. Examples of such assemblies include lipid rafts, and postsynaptic domains (PSDs) at excitatory and inhibitory synapses in neurons. PSDs are rich in scaffolding proteins that can transiently trap transmembrane neurotransmitter receptors, thus localizing them at specific spatial positions. Hence, PSDs play a key role in determining the strength of synaptic connections and their regulation during learning and memory. Recently, a two-dimensional (2D) diffusion-mediated aggregation model of PSD formation has been developed in which the spatial locations of the clusters are determined by a set of fixed anchoring sites. The system is kept out of equilibrium by the recycling of particles between the cell membrane and interior. This results in a stationary distribution consisting of multiple clusters, whose average size can be determined using an effective mean-field description of the particle concentration around each anchored cluster. In this paper, we derive corrections to the mean-field approximation by applying the theory of diffusion in singularly perturbed domains. The latter is a powerful analytical method for solving two-dimensional (2D) and three-dimensional (3D) diffusion problems in domains where small holes or perforations have been removed from the interior. Applications range from modeling intracellular diffusion, where interior holes could represent subcellular structures such as organelles or biological condensates, to tracking the spread of chemical pollutants or heat from localized sources. In this paper, we take the bounded domain to be the cell membrane and the holes to represent anchored clusters. The analysis proceeds by partitioning the membrane into a set of inner regions around each cluster, and an outer region where mean-field interactions occur. Asymptotically matching the inner and outer stationary solutions generates an asymptotic expansion of the particle concentration, which includes higher-order corrections to mean-field theory that depend on the positions of the clusters and the boundary of the domain. Motivated by a recent study of light-activated protein oligomerization in cells, we also develop the analogous theory for cluster formation in a three-dimensional (3D) domain. The details of the asymptotic analysis differ from the 2D case due to the contrasting singularity structure of 2D and 3D Green's functions.

The Influence of the Variable Wettability Characteristics of Layers on the Transport of Nanoparticles in the Context of Drug Delivery in Skin Structures.

The rapid development of nanotechnology has offered the possibility of creating nanosystems that can be used as drug carriers. The use of such carriers offers real opportunities for the development of non-invasive drug delivery through skin structures. However, in addition to the ability to create suitable nanocarriers, it is also necessary to know how they move through dermal layers. The human skin consists of layers with different wettability characteristics, which greatly complicates how introduced substances move through it. In this work, an experimental study of the diffusion process of nanoparticles through partitions with different wettability properties was carried out. Conventional diffusion tests using Franz chambers were used for this purpose. We quantified how the wettability of the barrier, the number of layers, and their mutual configuration affect the transport of nanoparticles. Based on the results, an analysis of the phenomena taking place, depending on the wettability of the partition, was carried out. A model relationship was also proposed to determine the effective diffusion coefficient, taking into account the influence of the wettability and porosity of the barrier.

In situ pulmonary mucus hydration assay using rotational and translational diffusion of gold nanorods with polarization-sensitive optical coherence tomography.

Significance: Assessing the nanostructure of polymer solutions and biofluids is broadly useful for understanding drug delivery and disease progression and for monitoring therapy. Aim: Our objective is to quantify bronchial mucus solids concentration (wt. %) during hypertonic saline (HTS) treatment in vitro via nanostructurally constrained diffusion of gold nanorods (GNRs) monitored by polarization-sensitive optical coherence tomography (PS-OCT). Approach: Using PS-OCT, we quantified GNR translational (DT) and rotational (DR) diffusion coefficients within polyethylene oxide solutions (0 to 3 wt. %) and human bronchial epithelial cell (hBEC) mucus (0 to 6.4 wt. %). Interpolation of DT and DR data is used to develop an assay to quantify mucus concentration. The assay is demonstrated on the mucus layer of an air-liquid interface hBEC culture during HTS treatment. Results: In polymer solutions and mucus, DT and DR monotonically decrease with increasing concentration. DR is more sensitive than DT to changes above 1.5 wt. % of mucus and exhibits less intrasample variability. Mucus on HTS-treated hBEC cultures exhibits dynamic mixing from cilia. A region of hard-packed mucus is revealed by DR measurements. Conclusions: The extended dynamic range afforded by simultaneous measurement of DT and DR of GNRs using PS-OCT enables resolving concentration of the bronchial mucus layer over a range from healthy to disease in depth and time during HTS treatment in vitro.

Travelling waves for a fast reaction limit of a discrete coagulation-fragmentation model with diffusion and proliferation.

We study traveling wave solutions for a reaction-diffusion model, introduced in the article Calvez et al. (Regime switching on the propagation speed of travelling waves of some size-structured myxobacteriapopulation models, 2023), describing the spread of the social bacterium Myxococcus xanthus. This model describes the spatial dynamics of two different cluster sizes: isolated bacteria and paired bacteria. Two isolated bacteria can coagulate to form a cluster of two bacteria and conversely, a pair of bacteria can fragment into two isolated bacteria. Coagulation and fragmentation are assumed to occur at a certain rate denoted by k. In this article we study theoretically the limit of fast coagulation fragmentation corresponding mathematically to the limit when the value of the parameter k tends to + ∞ . For this regime, we demonstrate the existence and uniqueness of a transition between pulled and pushed fronts for a certain critical ratio θ ⋆ between the diffusion coefficient of isolated bacteria and the diffusion coefficient of paired bacteria. When the ratio is below θ ⋆ , the critical front speed is constant and corresponds to the linear speed. Conversely, when the ratio is above the critical threshold, the critical spreading speed becomes strictly greater than the linear speed.

A psycholinguistic study of intergroup bias and its cultural propagation.

Intergroup bias is the tendency for people to inflate positive regard for their in-group and derogate the out-group. Across two online experiments (N = 922) this study revisits the methodological premises of research on language as a window into intergroup bias. Experiment 1 examined (i) whether the valence (positivity) of language production differs when communicating about an in- vs. out-group, and (ii) whether the extent of this bias is influenced by the positivity of input descriptors that were initially presented to participants as examples of how an in-group or out-group characterize themselves. Experiment 2 used the linear diffusion chain method to examine how biases are transmitted through cultural generations. Valence of verbal descriptions were quantified using ratings obtained from a large-scale psycholinguistic database. The findings from Experiment 1 indicated a bias towards employing positive language in describing the in-group (exhibiting in-group favoritism), particularly in cases where the input descriptors were negative. However, there was weak evidence for increased negativity aimed at the out-group (i.e., out-group derogation). The findings from Experiment 2 demonstrated that in-group positivity bias propagated across cultural generations at a higher rate than out-group derogation. The results shed light on the formation and cultural transmission of intergroup bias.

A simple method to prepare anion exchange membrane by PVA/EVOH/MIDA for acid recovery by diffusion dialysis.

Given the substantial environmental pollution from industrial expansion, environmental protection has become particularly important. Nowadays, anion exchange membranes (AEMs) are widely used in wastewater treatment. With the use of polyvinyl alcohol (PVA), ethylene-vinyl alcohol (EVOH) copolymer, and methyl iminodiacetic acid (MIDA), a series of cross-linked AEMs were successfully prepared using the solvent casting technique, and the network structure was formed in the membranes due to the cross-linking reaction between PVA/EVOH and MIDA. Fourier transform infrared spectrometer, X-ray photoelectron spectroscopy, scanning electron microscopy, and transmission electron microscopy were used to analyze the prepared membranes. At the same time, its comprehensive properties which include water uptake, linear expansion rate, ion exchange capacity, thermal stability, chemical stability, and mechanical stability were thoroughly researched. In addition, diffusion dialysis performance in practical applications was also studied in detail. The acid dialysis coefficient (UH+) ranged from 10.2 to 35.6 × 10-3 m/h. Separation factor (S) value ranged from 25 to 38, which were all larger than that of the commercial membrane DF-120 (UH+: 8.5 × 10-3 m/h, S: 18.5). The prepared membranes had potential application value in acid recovery.

Electrochemical Characterization of Diffusion in Polymeric vs. Monomeric Solvents.

Polymer electrolyte was used as a medium for testing the performance of microband electrodes under conditions of linear diffusion. Cyclic voltammetry (CV) and chronoamperometry (CA) experiments were performed in a highly viscous medium, where diffusion rates are much slower than in fluid solutions. The log i vs. log v (CV) or log i vs. log t (CA) relationships with the current equation confirmed the existence of such conditions, yielding slope values that were lower than the expected 0.5. This could indicate an impure linear diffusion profile, i.e., some contribution from radial diffusion (edge effects). However, the desired value of 0.5 was obtained when performing these tests in monomeric solvents of similar viscosities, such as glycerol or propylene glycol. These results led to the conclusion that the current equations, which are based on Fick's laws, may not be applicable for polymer electrolytes, where various obstructions to free diffusion result in a more complicated process than for monomeric solvents.

Coarse-grained modeling and dynamics tracking of nanoparticles diffusion in human gut mucus.

The gastrointestinal (GI) tract's mucus layer serves as a critical barrier and a mediator in drug nanoparticle delivery. The mucus layer's diverse molecular structures and spatial complexity complicates the mechanistic study of the diffusion dynamics of particulate materials. In response, we developed a bi-component coarse-grained mucus model, specifically tailored for the colorectal cancer environment, that contained the two most abundant glycoproteins in GI mucus: Muc2 and Muc5AC. This model demonstrated the effects of molecular composition and concentration on mucus pore size, a key determinant in the permeability of nanoparticles. Using this computational model, we investigated the diffusion rate of polyethylene glycol (PEG) coated nanoparticles, a widely used muco-penetrating nanoparticle. We validated our model with experimentally characterized mucus pore sizes and the diffusional coefficients of PEG-coated nanoparticles in the mucus collected from cultured human colorectal goblet cells. Machine learning fingerprints were then employed to provide a mechanistic understanding of nanoparticle diffusional behavior. We found that larger nanoparticles tended to be trapped in mucus over longer durations but exhibited more ballistic diffusion over shorter time spans. Through these discoveries, our model provides a promising platform to study pharmacokinetics in the GI mucus layer.

Spatially controlled diffusion range of tumor-associated angiogenic factors to develop a tumor model using a microfluidic resistive circuit.

Developing a tumor model with vessels has been a challenge in microfluidics. This difficulty is because cancer cells can overgrow in a co-culture system. The up-regulation of anti-angiogenic factors during the initial tumor development can hinder neovascularization. The standard method is to develop a quiescent vessel network before loading a tumor construct in an adjacent chamber, which simulates the interaction between a tumor and its surrounding vessels. Here, we present a new method that allows a vessel network and a tumor to develop simultaneously in two linked chambers. The physiological environment of these two chambers is controlled by a microfluidic resistive circuit using two symmetric long microchannels. Applying the resistive circuit, a diffusion-dominated environment with a small 2-D pressure gradient is created across the two chambers with velocity <10.9 nm s-1 and Péclet number <6.3 × 10-5. This 2-D pressure gradient creates a V-shaped velocity clamp to confine the tumor-associated angiogenic factors at pores between the two chambers, and it has two functions. At the early stage, vasculogenesis is stimulated to grow a vessel network in the vessel chamber with minimal influence from the tumor that is still developed in the adjacent chamber. At the post-tumor-development stage, the induced steep concentration gradient at pores mimics vessel-tumor interactions to stimulate angiogenesis to grow vessels toward the tumor. Applying this method, we demonstrate that vasculogenic vessels can grow first, followed by stimulating angiogenesis. Angiogenic vessels can grow into stroma tissue up to 1.3 mm long, and vessels can also grow into or wrap around a 625 µm tumor spheroid or a tumor tissue developed from a cell suspension. In summary, our study suggests that the interactions between a developing vasculature and a growing tumor must be controlled differently throughout the tissue development process, including at the early stage when vessels are still forming and at the later stage when the tumor needs to interact with the vessels.

Diffusion and Oligomerization States of the Muscarinic M1 Receptor in Live Cells─The Impact of Ligands and Membrane Disruptors.

G protein-coupled receptors (GPCRs) are a major gateway to cellular signaling, which respond to ligands binding at extracellular sites through allosteric conformational changes that modulate their interactions with G proteins and arrestins at intracellular sites. High-resolution structures in different ligand states, together with spectroscopic studies and molecular dynamics simulations, have revealed a rich conformational landscape of GPCRs. However, their supramolecular structure and spatiotemporal distribution is also thought to play a significant role in receptor activation and signaling bias within the native cell membrane environment. Here, we applied single-molecule fluorescence techniques, including single-particle tracking, single-molecule photobleaching, and fluorescence correlation spectroscopy, to characterize the diffusion and oligomerization behavior of the muscarinic M1 receptor (M1R) in live cells. Control samples included the monomeric protein CD86 and fixed cells, and experiments performed in the presence of different orthosteric M1R ligands and of several compounds known to change the fluidity and organization of the lipid bilayer. M1 receptors exhibit Brownian diffusion characterized by three diffusion constants: confined/immobile (∼0.01 µm2/s), slow (∼0.04 µm2/s), and fast (∼0.14 µm2/s), whose populations were found to be modulated by both orthosteric ligands and membrane disruptors. The lipid raft disruptor C6 ceramide led to significant changes for CD86, while the diffusion of M1R remained unchanged, indicating that M1 receptors do not partition in lipid rafts. The extent of receptor oligomerization was found to be promoted by increasing the level of expression and the binding of orthosteric ligands; in particular, the agonist carbachol elicited a large increase in the fraction of M1R oligomers. This study provides new insights into the balance between conformational and environmental factors that define the movement and oligomerization states of GPCRs in live cells under close-to-native conditions.

Direct visualization of pyrene diffusion in polyethylene and polyoxymethylene passive samplers.

While passive sampling of ultra-low aqueous concentrations of hydrophobic organic compounds in environmental aqueous media has emerged as a promising analytical technique, there is a lack of good understanding of the fundamental diffusive processes. In this research, we used a fluorophore, pyrene, as a model compound to track diffusion in polymers through absorption and environmental media exchange processes. We directly tracked the penetration of pyrene into polyethylene (PE) and polyoxymethylene (POM) rods during absorption from water by sectioning the rod after different stages of absorption and observing the fluorescence signal through a microscope. Diffusion profiles of pyrene in polymers were simulated by numerical integration of Fickian diffusion. The results indicated that the uptake process in PE is governed by Fick's law and the absorption and desorption kinetics are similar in this polymer. However, the observed uptake profiles of pyrene in POM were non-Fickian and the release kinetics out of POM was slower compared to uptake into the polymer. We show that slower desorption from POM makes corrections for nonequilibrium using performance reference compounds (PRCs) problematic for deployments in water or sediment where there is significant advection. However, for static sediment deployments, the overall kinetics of exchange is controlled by slow transport through sediment and the hysteretic behavior of POM may not preclude the use of PRCs to interpret equilibrium status.

Radiopharmaceutical transport in solid tumors via a 3-dimensional image-based spatiotemporal model.

Lutetium-177 prostate-specific membrane antigen (177Lu-PSMA)-targeted radiopharmaceutical therapy is a clinically approved treatment for patients with metastatic castration-resistant prostate cancer (mCRPC). Even though common practice reluctantly follows "one size fits all" approach, medical community believes there is significant room for deeper understanding and personalization of radiopharmaceutical therapies. To pursue this aim, we present a 3-dimensional spatiotemporal radiopharmaceutical delivery model based on clinical imaging data to simulate pharmacokinetic of 177Lu-PSMA within the prostate tumors. The model includes interstitial flow, radiopharmaceutical transport in tissues, receptor cycles, association/dissociation with ligands, synthesis of PSMA receptors, receptor recycling, internalization of radiopharmaceuticals, and degradation of receptors and drugs. The model was studied for a range of values for injection amount (100-1000 nmol), receptor density (10-500 nmol•l-1), and recycling rate of receptors (10-4 to 10-1 min-1). Furthermore, injection type, different convection-diffusion-reaction mechanisms, characteristic time scales, and length scales are discussed. The study found that increasing receptor density, ligand amount, and labeled ligands improved radiopharmaceutical uptake in the tumor. A high receptor recycling rate (0.1 min-1) increased radiopharmaceutical concentration by promoting repeated binding to tumor cell receptors. Continuous infusion results in higher radiopharmaceutical concentrations within tumors compared to bolus administration. These insights are crucial for advancing targeted therapy for prostate cancer by understanding the mechanism of radiopharmaceutical distribution in tumors. Furthermore, measures of characteristic length and advection time scale were computed. The presented spatiotemporal tumor transport model can analyze different physiological parameters affecting 177Lu-PSMA delivery.

MAM-E: Mammographic Synthetic Image Generation with Diffusion Models.

Generative models are used as an alternative data augmentation technique to alleviate the data scarcity problem faced in the medical imaging field. Diffusion models have gathered special attention due to their innovative generation approach, the high quality of the generated images, and their relatively less complex training process compared with Generative Adversarial Networks. Still, the implementation of such models in the medical domain remains at an early stage. In this work, we propose exploring the use of diffusion models for the generation of high-quality, full-field digital mammograms using state-of-the-art conditional diffusion pipelines. Additionally, we propose using stable diffusion models for the inpainting of synthetic mass-like lesions on healthy mammograms. We introduce MAM-E, a pipeline of generative models for high-quality mammography synthesis controlled by a text prompt and capable of generating synthetic mass-like lesions on specific regions of the breast. Finally, we provide quantitative and qualitative assessment of the generated images and easy-to-use graphical user interfaces for mammography synthesis.

A Study of the Drift Phenomena of Gate-Functionalized Biosensors and Dual-Gate-Functionalized Biosensors in Human Serum.

In this paper, we study the drift behavior of organic electrochemical transistor (OECT) biosensors in a phosphate-buffered saline (PBS) buffer solution and human serum. Theoretical and experimental methods are illustrated in this paper to understand the origin of the drift phenomenon and the mechanism of ion diffusion in the sensing layer. The drift phenomenon is explained using a first-order kinetic model of ion adsorption into the gate material and shows very good agreement with experimental data on drift in OECTs. We show that the temporal current drift can be largely mitigated using a dual-gate OECT architecture and that dual-gate-based biosensors can increase the accuracy and sensitivity of immuno-biosensors compared to a standard single-gate design. Specific binding can be detected at a relatively low limit of detection, even in human serum.

Analytical solution of fuzzy heat problem in two-dimensional case under Caputo-type fractional derivative.

This work aims to investigate the analytical solution of a two-dimensional fuzzy fractional-ordered heat equation that includes an external diffusion source factor. We develop the Sawi homotopy perturbation transform scheme (SHPTS) by merging the Sawi transform and the homotopy perturbation scheme. The fractional derivatives are examined in Caputo sense. The novelty and innovation of this study originate from the fact that this technique has never been tested for two-dimensional fuzzy fractional ordered heat problems. We presented two distinguished examples to validate our scheme, and the solutions are in fuzzy form. We also exhibit contour and surface plots for the lower and upper bound solutions of two-dimensional fuzzy fractional-ordered heat problems. The results show that this approach works quite well for resolving fuzzy fractional situations.

Modelling non-local cell-cell adhesion: a multiscale approach.

Cell-cell adhesion plays a vital role in the development and maintenance of multicellular organisms. One of its functions is regulation of cell migration, such as occurs, e.g. during embryogenesis or in cancer. In this work, we develop a versatile multiscale approach to modelling a moving self-adhesive cell population that combines a careful microscopic description of a deterministic adhesion-driven motion component with an efficient mesoscopic representation of a stochastic velocity-jump process. This approach gives rise to mesoscopic models in the form of kinetic transport equations featuring multiple non-localities. Subsequent parabolic and hyperbolic scalings produce general classes of equations with non-local adhesion and myopic diffusion, a special case being the classical macroscopic model proposed in Armstrong et al. (J Theoret Biol 243(1): 98-113, 2006). Our simulations show how the combination of the two motion effects can unfold. Cell-cell adhesion relies on the subcellular cell adhesion molecule binding. Our approach lends itself conveniently to capturing this microscopic effect. On the macroscale, this results in an additional non-linear integral equation of a novel type that is coupled to the cell density equation.

Modeling of brain efflux: Constraints of brain surfaces.

Fluid efflux from the brain plays an important role in solute waste clearance. Current experimental approaches provide little spatial information, and data collection is limited due to short duration or low frequency of sampling. One approach shows tracer efflux to be independent of molecular size, indicating bulk flow, yet also decelerating like simple membrane diffusion. In an apparent contradiction to this report, other studies point to tracer efflux acceleration. We here develop a one-dimensional advection-diffusion model to gain insight into brain efflux principles. The model is characterized by nine physiological constants and three efflux parameters for which we quantify prior uncertainty. Using Bayes' rule and the two efflux studies, we validate the model and calculate data-informed parameter distributions. The apparent contradictions in the efflux studies are resolved by brain surface boundaries being bottlenecks for efflux. To critically test the model, a custom MRI efflux assay measuring solute dispersion in tissue and release to cerebrospinal fluid was employed. The model passed the test with tissue bulk flow velocities in the range 60 to 190 [Formula: see text]m/h. Dimensional analysis identified three principal determinants of efflux, highlighting brain surfaces as a restricting factor for metabolite solute clearance.

A mathematical model of tissue axial and radial diffusion in the microvasculature for intravascular microscopy and phosphorescence quenching data.

This study aims to extend earlier Krogh Cylinder Models of an oxygen profile by considering axial diffusion and analytically solving Fick's Law Partial Differential Equation with novel boundary conditions via the separation of variables. We next prospectively collected a total of 20 animals, which were randomly assigned to receive either fresh or two-week-old stored red blood cell (RBC) transfusions and PQM oxygen data were measured acutely (90 min) or chronically (24 h). Transfusion effects were evaluated in vivo using intravital microscopy of the dorsal skinfold window chamber in Golden Syrian Hamsters. Hamsters were initially hemorrhaged by 50% of total blood volume and resuscitated 1-h post hemorrhage. PQM data were subsequently collected and fit the derived 2D Krogh cylinder model. Systemic hemodynamics (mean arterial pressure, heart rate) were similar in both pre and post-transfusion with either stored or fresh cells. Transfusion with stored cells was found to impair axial and radial oxygen gradients as quantified by our model and consistent with previous studies. Specifically, we observed a statistically significant decrease in the arteriolar tissue radial oxygen gradient after transfusion with stored RBCs at 24 h compared with fresh RBCs (0.33 ± 0.17 mmHg µ m-1 vs, 0.14 ± 0.12 mmHg µ m-1; p = 0.0280). We also observed a deficit in the arteriolar tissue oxygen gradient (0.03 ± 0.01 mmHg µ m-1 fresh vs. 0.018 ± 0.007 mmHg µ m-1 stored; p = 0.0185). We successfully derived and validated an analytical 2D Krogh cylinder model in an animal model of microhemodynamic oxygen diffusion aberration secondary to storage lesions.

A hard sphere model for single-file water transport across biological membranes.

We use Gürsey's statistical mechanics of a one-dimensional fluid to find a formula for the P f / P d ratio in the transport of hard spheres across a membrane through a narrow channel that can accommodate molecular movement only in single file. P f is the membrane permeability for osmotic flow and P d the permeability for exchange across the membrane in the absence of osmotic flow. The deviation of the ratio from unity indicates the degree of cooperative transport relative to ordinary diffusion of independent molecules. In contrast to an early idea that P f / P d must be equal to the number of molecules in the channel, regardless of the physical nature of the interactions among the molecules, we find a functional dependence on the fractional occupancy of the length of the channel by the hard spheres. We also attempt a random walk calculation for P d individually, which gives a result for P f as well when combined with the ratio.